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IITRODUOTION 



Previous Investigation devoted to the study of the 
turlsulent "boundary layer and drag of wings and figures of 
revolution were "based essentially on the assumption of the 
power law of velocity distrihution in the "boundary layer, 
the exponent n in the velocity distrihution formula 

u = (IJ (1) 

(where u is the velocity in the "boundary layer, Ug = 
the velocity at the outer limit of the 'boundary layer, 
y = the ordinate taken normal to the wall, and 6 = the 
thickness of the "boundars' layer) "being taken equal to l/7. 
Tising this method the prohlem of the cylindrical wing was 
solved "b^ Miller (reference l) and that of the airship 
"body "by Millikan (reference 2), 

However, as various tunnel tests have shown (see 
Nikuradse, reference 2), the form of curve (l) varies with 
the Reynolds Number in such a manner that to a first ap** 
proximation the exponent successively assumes the values 
l/s, 1/9, 1/10, etc., as the Eeynolds Humher increases* 

A certain forward step in investigating the effect 
of a change in n • with Reynolds Numher on the "boundary 
layer and the drag of figures of revolution was taken in 
1932 hy K. K. Pedlaevsky (reference 14) who, making use of 
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the test data of Nikuradse for smooth tuhes, found an ap- 
proxima.te relation hetween the exponent n and the EeynoXds 
Numher formed from the Telocity at infinity and the length 
of the body (Re = VL/u)* 

The aboYe-mentioned method possessed, however , two 
unavoidable theoretical defects, namely: l) the computa-* 
tion, on the basis of Nikuradse tests, of the relation 
n = f(Ee) was very approximate and 2) the Tariation of 
n along the body at constant Ee , was neglected. The 
first defect may only be removed by conducting numerous 
tedious wind-tunnel tests with figures of revolution, sim- 
ilar to the tests of Hikuradse, and these would hardly 
give satisfactory results on account of the smallness of 
the Eeynolds Humbers attainable . in our wind tunnels. With 
rego.rd to taking account of the variation of n along the 
length of the body due to the gradual increase in the lo-* 
cal Eeynolds Number along the body, we should be confront- 
ed with very great mathematical difficulties in the solu- 
tion of the boundary-layer equation in attempting to find 
and to substitute in it some relation between n and the 

6 u 

Eeynolds lumber Eg = 

With the above considerations in mind, we proceeded 
to work out a method for applying to the solution of the 
problem of the turbulent boundary layer and drag of stream- 
liixe bodies, the universal logarithmic law of velocity 
distribution theorotically developed by Prandtl (reference 
5) and K^rman (reference 6) in 1929-1930. As will appear 
below, the universal logarithmic law by the very prin- 
ciple on which it is based avoid.s both the difficulties 
mentioned above and offers, if not a final solution, at 
least one that much more closely approaches actual condi-^ 
tions» 

Tile universal logarithmic formula of Zarman was. used, 
as giving better agreement with experiment and as being 
derived from more rigorous assumptions than the correspond- 
ing formula of Hikuradso-Prandtl ♦ As will be shown below, 
however, our equations are entirely applicable also to the 
case where the Uikuradse-Prandtl velocity distribution for- 
mula is taken as a starting point since the difference 
will evidence itself only in the change of certain con- 
stant coefficients . 

In the present paper we shall consider a figure of 
revolution, so that the formulas applicable to the more 
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sinple cases as, for example, a wing or a flat plate will 
follow from our equatioiis as corollaries, 

Por checking the results of our theory, we made use 
of the data derived fron the excellent tests of Freeman on 
a l/40-scale raodel of the airship "Akron" conducted in the 
large N«A«G»A, wind tunnel • 

At the time when our work was already completed in 
its essentials, we ohtained a preliminary report on the 
work of Moore (reference 7), in whicii the prohlem of a 
figure of revolution i s' so Ived hy application of the 
Nikuradse-»Prandtl law* Our equations, however, appear .to 
he more general in the sense that they are applicahre to 
various "bodies and appear moreover to he exact solutions 
of the integral relation of Karman which lies at the has is 
of the houndary-layer theory, since we are the first - so 
far as we know - to take full account of the curvature 
of the concentric elements of the houndary layer, i»e», 
their va,rious distances from the axis of the hcdy of rev*^ 
olutlon and .to show numerically tha.t the variation is 
small and does not appreciaoly affect the final result - ■ 
thus rigorously proving that the curvature nay he' neg-^ 
iected. In addition we give, in our opinion, a more ra- 
tional procedure for the numerics&i integration of the 
equations ohtained than does Moore. The greatest differ- 
ence, however, "between our work and the similar work of 
other authors lies in hrlnging to light the large devia^ 
tions of the results of application of the Karman theory 
in its origina.1 form from the results of experiments* In 
the second part of this paper we have attempted to explain 
these discrepancies hy suhjecting .the logarithmic theory 
to a general criticism. We have thus succeeded in deriv- 
ing theoretically a sufficiently justifiahle method of 
rendering the logarithmic theory more precise, after 
which it was found to he in excellent agreement with ex— 
periment^ 

In the first part we shall derive the fundamental 
equation for a "body of revolution according to the K^rm^n 
theory in its original form, and in the second part we 
shall give all the comparisons of the results of tests 
with the modified theory^ 
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PAR T I 

1. OTIYSESAL LOaARITHMIG FORMULAS 

Before proceeding with the prohlem itself, we shall 
hriefly review the essentials of the Earman. theory* 

Earman starts out from the following two fundamental 
assunptions: l) Both the turbulent raixing characteristics 
(the mean nixing length) and the turbulent friction at a 
given point are completely deterniined hy the characterise* 
tic nagnitudes of the mean flow at the gi-ren point - in 
particular, "by the deriYatiVes of the mean flow; 2) The 
turbulent flow pulsation patterns are similar:, i.e., dif- 
fer only in the time and distance scales 8.nd do not de- 
pend on the Yiscosity. 

. As the tests of Stanton and the later ones of Fritch 
have shown, the velocity-distribution patterns at pipe sea- 
tions having different wall roughness - i.e., having a dif-* 
ferent shear stress at the walls - appear to be entirely 
identical provided the difference between the maximum ve-* 
locity at the channel axis and the velocity at some other 

/tT 

point is expressed in terms of J'p' ^^^^^re is the 

shearing stress or skin friction per unit area at the wall 
and p is the mass density: 




In the above z is the distance from the pipe axis at 
right angles to the mean flow direction and r is the 
pipe radius. Thus the expression at the left is a func--* 
tion only of the poordin8.tes and does not depend on the 
viscosity. This, naturally, is true only for regions not 
too n e a r the wall. 

The above suggests the existence of actual mechanical 
similitude between turbulent flows having different dynamo- 
teal boundary conditions (shear at the wall). 

If the form of the function f(z/r) in expression 
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i^) is foxind, the prol3lem is easily solred with regard to 
the velocity distrihution iii the pipe, the maximum Teloc- 
ity *^inax' pips resistance. 

The most important feature of formula (3) is that it 
is true for any sufficiently large Eeynolds lumber; i»e», 
that the form of the function f(z/r) does not depend on 
the Reynolds Humher. For small Reynolds Humhers, where 
the effect of the riscosity shows up more strongly, formu*^ 
la (3) has less Yalidity, 

The prohlem of finding the function f(z/r) was 
solved hy Kd^rm^n (reference 4) "by applying the ahove as-^ 
sumptions. of mechanical similitude of the flow pulsation 
patterns and the absence of the effect of viscosity to the 
fundamental hydrodynamie equations (reduced to the Helm- 
holtz form on the assumption of the existence of a flow 
function for the actual velocities):. He ohtalned a new 
formula for the mixing length I : 

u " 

where and u^ are derivatives of the velocity with 

respect .to the coordinate at right angles to the mean flow 
and K is a universal constant. Eliminating I "by the 
known expression 

T ^ p u»^ (4) 

'(\?hich, as Karm.an has shown, also appears as the result of 
the assuiiiption of similitude of the pulsation field) where 
T is the shearing stress at any point of the flow, and 
assumin^^ T to vary linearly with the distance z from 
the axis of the pipe; that is, = ^ A Karman obtains 

the fundamental .differential equation 

u » ^ ■ z 



. - (5) 
Ti O r 



Upon integrating, with the corresponding choice of 
the constants of integration (at z ^ 0 u = Uj^^g^^c 

2 =: r u * " Ts) , there is obtained ; 

"^Wis expression is a result of the requirement of dynamic 
eq^uilibrium of the concentric walls of the fluid cylinders 
for the established mean fl»w. 
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i,e.^ the form of the function f (z/r) is completely de-* 
teriiiined. 

Introducing, in place of the distance z from, the' 
axis of the pipe, the distance y from the wall (z - 
r - y) from (6), we find for the ratio ii/*^niax 



= 1 -f 



max 



1 

K u" 



max 



As . nay he seen, the ratio ^l^rs 



is not suhject to the 



power law* It is interesting to note that as the Reynolds 
Mumher increases - that is, with decreasing value of ^'^qIp, 
curve (o 0 changes its form, the curves ohtained closely 
approaching the power law with consecutively decreasing 
exponents - i/7» l/S, l/9f etc. 

The dependence of the velocity profile on the Reynolds 
Humher is thus already contained in formula (6'), It is 
not difficult to show that etjuatlon (6') has a discontinui- 
ty at y = r; that is, on the axis of the pipe - a result 
which is due to the fact that derivatives of u with re-* 
spect to y higher than the second were neglected. Ihe 
effect of this is only slight since the break in the smooth- 
ness of the curve (6) at y = r is very small. At y ~ 0, 
that is, at the wall itself, the function (6) gives infi- 
nitely large velocities u of logarithmic order. This is 
also explainahle theoretically since the entire theory 
ceases to he true near the wall as a result of the strong 
increase in the effect of the viscosity. The latter re- 
sults in the formation of a thin laminar layer at the wall, 
the existence of such a layer having heen verified hy nu- 
merous tests* The thickness of this layer and the veloci- 
ty just outside it can only depend on the friction at the 
wall, the viscosity, and the fluid density. On the hasis 



of dimensional considerations 
pressed as 



this thickness may he ex- 



*lam 



= 3 



(7) 
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and the velocity "^S^gj^yj^ J'^-st outside the layer as 



Sy the use of theahove expressions Karman was ahle 
to derive the law for the resistance of a round pipe« 
ABSXimlng law (6) to he valid from the outer houndary of 
the laniiiar layer to tho center of the pipe, that is, from 
y ^ ^\B,jr^ y - 37, the value of the velocity at the 

pipe axis may he found from expression (6), written for 

y = Slam- 



max - ^^6 



nam K p L " 3: / ^ r 

Since ^i^m very small hy comparison with r, 

this expression may he simplified hy extracting the approx- 
im.ate root 



^ ^ 1 am ^ n ^ lam 



r " Si- 

There is then ohtained 

1 /tT r 2r 



- U;i„ -J- -I Mn 1 

1 am 



(9) 



Substituting expressions (?) and (8), we obtain; 



"^^max " K 




+ c/ (10) 



where 0 includes the constants A, B, etc* 

Substituting in expression (9) the Reynolds iNTumber of 
the pipe R = and the nondimeiisional frlctional 

T*o 

coefficient Gj> = ^ obtain: 

5 P ^ max 



^ In (EycTT + d (11) 
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where G * == G + In , 

This dependence of the resistance coefficient on the 
Reynolds ITumher agrees excellently urith tests oyer a wide 
range of Edynolds lumlDers. 

In 1930, Prandtl ohtained a -aniversal logarithiniG law 
of velocity distrihution somewhat different from the law 
discussed ahove , of Karman. Prandtl considered the infi- 
nite flow about a fixed wall and assumed the mixing length 
to "be proportional to the distance y from the wall 
i»e*, in place of expression (4)* he had 

I - E y (4» ) 



This gave, after integrating the expression t= p I 



u 



= - 1 11 y + c 



Hikuradse, on. the tasis of his tests, ohtained the 
formula 

" ^aax ^ ^ i In 1 (6") 




which may easily "be, ohtained from the Prandtl formula hy 
applying it to a pipe, excluding the constant c, and 
using the condition u = u^g^^ at y = r^ This is the 
formula also used "by Moore* 

It is interesting to note that if expression (6") is 
applied to the laminar layer at the wall as was done ahove 
with formula (6), we ohtain the same expression (lO) hut 
with the constanx G differently comhined out of the un- 
known coefficients A, B, etc. Thus the same resistance 
law is ohtainedo The constants 0' and K may he de-» 
termined only experimentally. Up to the present these 
constants have hoen dotermined by Karman and likuradse en- 
tirely "by experiments on the dependence of the resistance 
of pipes on the Reynolds Numher R, that is, hy compar-^ 
ing formula (ll) with tests. In doing this It was assumed 
that the universal constant E thus determined is iden- 
tical with the constant K that enters into the velocity 
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formulas (6) or (6"), In the second part of our work we 
shall .giye special consideration to this prohlem and show 
that actually this is far from "being true. For the pres- 
ent, we shall pass on to the derivation of the equation- 
for the solution of the hou-ndary^layer prohlem and the 
drag of a figure of revolution, using as a hasis the ve^ 

• locity and drag-^distrlhution formulas of . Karman in the 

■ general form in which he derived them. 

2, DEaiTATIOSr OF THE FJNDAMINTAL EQTJATIOl 



. In the previous section we considered the laws of 
flow and the resistance of a circular pipe as derived by 
Karman from the theory of mechanical similitude* The solu- 
tion of Karman refers to the so-called "internal prohlem." 
Our ol^ject in this work is to oh tain a solution of the 
"external prohlem," i.e., the case where in place of a 
flow enclosed within a cylindrical pipe we have the flow 
of an infinite mass of fluid ahout a "body. This hody in 
our ease will he a figure of revolution - for example, 
the hull of an airship. 2) fg shall not distinguish the case 
of motion of a hody in an undisturhed- medium (airship 
flight) from that of a turbulent stream flowing up to- a 
body fixed in position (aerodynamic tunnel) although these 
cases are not equivalent from the point of view of drag. 
The drag, as we know, depends chiefly on the character of 
the bqundary layer which in turn is determined exclusively 
by the Reynolds Numbor Vl/p computed from the velocity 
of the motion (or the velocity of the mean flow at infin--* 
ity) and some linear dimension of the body. The Eeynolds 
Number affects both the absolute thickness of the boundary 
layer and the position of transition of -the layer from the 
laminar to the turbulent condition* Since both of these 
factors also depend on the degree of turbulence of the 
flow in th-e tunnel, this degree o-f turbulence may be taken 
into account, in the final computation, by a correspond-^ 
ing increase in the Eeynolds Number. It is true that 
since up to the present we possess no method for the quan- 
titative determination of the degree of turbulence, we 
cannot as yet perform such a numerical operation. For 
this reason, strictly speaking, we may only study the mo-* 
tion of body in an undisturbed medium. This is the case 
of greatest practical interest. 



In this v/ork we shall consider only bodies with smooth 
surface . 
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We shall apply the conclusions - of the preyious 'section 
in their entirety to the external prohlem. It is therefore 
^necessary to consider the justif lability of this procedure. 

The turhulent houndary layer ahout a body situated in 
a flow is that region- in which the action of the turbulent 
friction destroys the Telocity distribution given by the 
potential flow. We haTe the same condition also in a cir- 
cular pipe. In this case these f rictional-^ f or ces give the 
velocity profile of expression (6^)* Since the latter was 
derived on tha assumption that the effect of the turbulent 
friction shows up throughout the pipe cross section, it 
follows that in this case the entire region within the pipe 
is a "rolled up" boundary layer whose thickness. is equal 
to the pipe radius and whose outer 'limit has become a sin-- 
gle straight line, namely, the axis of the pipe, The ve- 
locity at- this axis is, according to formula (10), a func- 
tion of the frictional intensity at the wall, i.e..,- a func- 
tion of the state of the boundary >layer . In the case of 
the externevl problem this velocity just outside the bound- 
ary layer will appea.r as the boundary condition determined 
by .the potential' flow outside the boundary layer. This is 
•where the difference lies between the internal and exter- 
nal problems. The velocity just out-side the boundary 
layer wo shall in what follows denote by Ug • 

Strictly speaking, we are not quite justified in pass- 
ing directly from the internal flow along the rectilinear 
wall of the pipe to the external flow a,long the curved - 
boundary of a wing or figure of revolution. It has be^n 
established by tests that the flow, for example, in diverg- 
ing and converging channels is not subject to the loga- 
rithmic law. It would be inore rigorous if the transition 
to the externa.l problem were made from the latter two 
cases. Since, however., we have up to the present, no so- 
lution of the problem of turbulent flow in diverging and 
converging channels, we are required to make use of the 
generally adopted method of transition to the external 
problem from straight pipes. 

Thus, denoting by 8 the thickness of the boundary 
layer, we shall have for the external case as the velocity 
distribution in the boundary layer 
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where y is takea normal to the surface of the body, 

jormula (11) giving the dependence of the frictional 
coefficient on the Eeynolds Humher, "becoines: 



:= In (Hg /Gj) + G ' (13) 



where the coefficient is the ratio of the local fric- 

tional intensity to the dynamic velocity just outside the 
"boundary layer *^ 

(14) 




and the Eeynolds Wumher is Gomputod from the velocity just 
outside the "boundary layer and the thickness of the lajrcr: 

u« 8 

As may he seen, in the case of the external prohlcm 
we "bring into consideration, in place of the knoTrn radius 
of the pipe the unknown thickness "of the houndary layer 
which, as is shown hy "both theory and experiment, varies 
along the hull« Thus, for a complete solution of the proh- 
lem, we must have still another condition that connects 
the variahle thickness 6 of the houndary layer with the 
other variahles, of which Ug is considered as given "by 

the assumption of a potential flow just outside the hOund-' 
ary layer. 

An additional relation of this kind appears to "be the 
well-known Integral rela,tion of Karman which we shall now 
considor. This relation, derived for two-dimensional, 
steady flow assumes, for a figure of revolution, the fol*^ 
lowing f orm : 



es J 2'^^' 'P'^ dy-ug J 2Trr'pudya - — J a-rrr 'dy-2Trr Tp 
0 'o 'o (15) 

^3?he local frictional intensity as rre shall see farther on, 

■crill vary along the hull of the body,. . , , 
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where s is the length along the hull of the figure of 
reTolut ion , measured from the nose. 

r', radius of a concentric element of the surface. 
(See fig. 1/) 

p, density of the fluid. 

y, distance from the surface measured ^ilong the 
normalo 

8, thickness of the houndary layer. 

U and Ug , mea.n velocities within and just outside 
the layer . 

p> static pressure at "boundary layer. 

r, radius of circular section of figure of rero-* 
lut ion, 

Tq , local shearing stress (or skin friction per 

unit of area) and the variahles have the follow*^ 
Ing functional relations anong each other : 

r' = r» (s,y); 6 == 6(s); u = u(s,y) 

^ Ug(s); p = p(s); r = r(s) 

Expression (16) thus gives the required additional 
relation expressing either as a. function of 6 a^nd s, 

( 6 , s ) , or 6 as a f unc t ion of and s , 6 ( , s ) , and 

the problem thus "becomes determinate. In equation (16) the 
left*^hand side represents the increase per second of the 
momentum of the fluid at an element of houndary layer of 
length ds and the right-^Iiand side gives the sum of the 
forces acting on the "boundary layer. This relation will 
he pro ve d h e 1 o w « 

Figure 1 shows a portion of a section of a figure of 
revolution (generating line ah) irith the boundary layer of 
thickness 6* Through point A on the body, pass a, conica^l 
surface normal to the surface of the body. Let the gener- 
ating arc ab rcce i ve an increment as = AA ' and through 
point A^ pass a conical surface A ^3^ normal to the sur^ 
face of the body. Through the circle of intersection of 
surfa.ce A'B* vilth the outer surface of the boundary layer 
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ve pass the "surf ace . of flow" a'i)V. We thus ohtain an 
annular element of a "tutie of flow" "bounded "by the surface 
of flO'T7 and the "body. To this element we may apply the 
" no me n turn 1 aw . " 

Let us denote the mass of fluid flowing per second 
through section AB of the "boundary layer "by Q, and the 
monentuii associated with this mass hy J. These nagnitudes 
are glTen hy 



Q = / 2Trr»pudy 



J - / 27Tr *pu dy 
o 

In passing through the element 
i t s nonentuin J increase "by 



ds , 

and 



the mass 



ds.. 



Q an d 
respec- 



d s 

d s d s 

tively; that is, the mass and the momentum of the fluid 
flowing through section A^B' will he 



J + 



I- ^= 



The increase in the nass —-^ ds can only he obtained hy 

9 s 

the entrance of a nass of the potential flow fluid through 
portion B Bi of the section of t he tuh e of flow^ Since 

the Telocity of this inflow is Ug the momentum of the 

fluid introduced through the entire section ABi will he 



J 4- 



^6 ? 



ds 



ds 



Our suri;a,ce of flow was so drawn that the section A VB ^ 
of the houndary layer is at the sane tine a section of a 
tuhe of flow» Thus the increase in the nomentun of the 
fluid through an ele:aent of the tuhe of flow will "be 



r.6 



J ^ M a, . J . If ds . ^ 2Trr^pudy 

o 

r /v5 



ds 



" "".^ fs Ij 2Tir'pudy 

0 



ds (a ) 
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This increase of nomentuin per second should, accord- 
ing to the noirientun law, he equal to the sun of the forces 
acting on the element of the tuhe of flow in the direction 
of the nonentun vector. If the static pressure at the left 

of the elenent is td and that at the right p + 12 ds the 

3 s 

force ^ eYidenced hy the drop in pressure, will, be 



dsy / SiTr ^dy-p / 2tTr'dy 



6 + 1^ ds 

OS 



J 2Trr »dyl ds 



The ninus sign is chosen since a positive pressure gradX'- 

ent rr— , corresponds to a force in the optiosite direction 
as 

to the no men turn vector. 

Breaking up the integral thus ohtained into two, we 

have 

6+ 1^ ds 

8 3S g 

- |2 / 2Trr »dy+?£ / 2Trr »dy 1 ds ^ - [ |£ f 2Trr'dylds 
t as j as J J ^°^.J j 

The second integral vanishes since it is infinitely small 
in comparison with the. first,: The resultant of the turhu- 
lent frictional forces on the element ds will .he 

- 2tt r Tq ds (c) 

Equating expression (a) to the sun of expressions 
(h) and (c) and dividing through hy ds , we obtain the 
integral relation (16) of Karman. 
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We liave assumed throu^fnout that the static pressure 
does not vary along the thickness of the houndary layers 
Aetuaily there exists a certain, thotigh very snail, pres-^ 
sure gradient in the "boundary layer (see, for example, 
reference 12), of which no quantitative account can as yet 
"be taken. Attempts to determine the normal forces re--* 
suiting from the pressure gradient normal to the surface 
hy equating them to the integral of the centrifugal forces 
of the particles of fluid moving on a given curvilinear 
contour in the houndary layer gave negative results the 
pressure gradient ohtained as a result of the centrifugal 
forces making up, for example, only half the value that is 
actually oliserved. 

Our problem now is to "bring equation (16) into rela-^ 
tion with expressions (iS) and (14). We shall first tran.s-^ 
form expression (15). The radius r' of the concentric 
elerient of the houndary layer, according to figure 1, is 
e qual t o 

r * = r + y cos G 

where G is the angle "between the tangent to the generat*^ 
Ing line of the "body of revolution and its axis- The 
quantity y cos G is generally neglected hy comparison 
with r; that is, r^ is considered equal to r and the 
effect of this simplification on the final result is not 
analysed. We shall not introduce this simplification in 
order that we may have the possibility of determining the 
effect of neglecting y cos 6 on the final result. ' Sub*^ 
stituting in expression (16) the value of r^, dividing 
by 2tt P and transferring to the left-hand side the mem- 
bers containing cos 6, we obtain: 

6 8 6 8 

I ru^dy-^U5 j rudy + / cos 6 u ydy^Ug ~ y^s 6u:ydy== 

0 6 ' 0 0 

6 

1 rST) P, V "^O 

= - i / (r+y COS e ) dy-r -^r- 
p d s J • p 

O 

Taking r and cos G outside the integral signs as being 
independent of y and performing all the possible differ-* 
entiations and inte gr a t i on s , we have : 
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9 P^a 3r /'^a B Br 

r 3- / u dy + / u dy - rug ^- / udy - Ug / udy + 

0 0 '0 p 

+ cos 0 I- / ,T#ydy+ i-Hs„§. / ydy-ug cps 9 |- •/ 



8 cos 6 r ^ 1 dp / ^ 1 T5 . . 

^ U5 — / uydy = p of (^^^"^ 2 a cos 9^ - r ~ (17) 



To the region outside the ■boundary layer we may apply 
the Bernoulli theorem and thus express the pressure gra-- 
dient entering into the expression on the right-hand side 
in terms of the velocity just outside the ■boundary layer. 
We have : 

P i^A^ . j. 

t> + : - constant 

2 

Differentiating with respect to s, we ohtain: 

1 dp 3ug 
" p " 3s" 

Grouping the terms in expression (17) and suhst i tuting the 
value of ^ --^^ 5 we ha.ve : 




Making use of the velocity distrihution law (12), we must 
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now compute the Integrals in expression (18). 



We shall, denote the expro^sioii 



V p 



by 



This 



expression, whioh appears as the "dynamic scale" of veloc-* 
ities in equation (12) and has the dimensions of a veloc- 
ity, we shall denote - following the suggestion of Profes- 
sor !• &• Loitsansky - as the "dynamical Yelocity." We 
shall further introduce a new rariahle Xi 



X = 1 



(19) 



noting that a variation in y from 0 to 6 corresponds 

to £1 variation in x from 0 to !• We. then ohtain from 
expression (12) the velocity and its square: 



U == Ub 



K 



[in X + 1 x] 



u^ = ui + 2 [In X X -^ 1] + -^s- (In X ^ X + l)' 



(20) 



We may no-cr compute all the integrals, "bearing in mind that 
in a^^^reemeiit Trith the change in variables (19)2 

y = 6 - 6 (1 - x)^ 

dy = 2S (1 - x) dx 

Without going through all the computations, we shall at once 
v/rite doiYn the final result: 



udy - 



J 



>8 



6 



\ • 6 E 9 irJ 



uydy = 



/ 23 V* 



) 



/'^8 ^ 1 .2 / 46 UgV* 601 



(21) 
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In the process of integration we meet with integrals 
with rest)ect to x, and it is useful to write these down: 



1 ) / In xdx = fx In x x 1 ^ = - 1 
. / 0 



2 ) / X In xdx = i ™— 



-,1 



1 
4 



3) 



2 ^ i In X 



In xdx = 



a 



"1 



4 ) / x*^ In xdx = I 



X 



4 1 



16 J, 



X 
9 



IB 



5) / In^xdx = [x In^x - 2x In x + 2x1^ = 2 
J 

0 



6) / X in^xdx = i^^i^--^ - 



la 



4 



1 
4 



7) / x^ lii^ xdx 



x^ In^ X 



2x^ In X ^ 2x^ 



t1 



27 J 



2_ 
27- 



8) 



In^ xdx =: I 



x^ In^ X In X 



4 



8 



x^ ■ 
32 



1^ 



All the s.hove integx'als are readily integrated hy parts 
noting that the lintt of the expressions In^ x, ^vhere 

n and m are positive integers, is sero as x approaches 
sero, as naj readily "be proved "by repeated application of 
the L VHo sr. i t al rul e . 



Before suDstitut ing the integrals (21) in expression 
(18) we shall eliiiiina.te from them 5. Using the • expressions 
for the friction coefficient 0^ and the Eeynolds Numher 
E.5 as given hy (14) a,nd (15), we ohtain from (13): 
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or 



V pa 



aiid ' introducing the . constant C * under the logeri thm sign, 
i.e., suhst i tut ing e^a = Cs, we ohtain on solving the 
above equation for 8 : 

KUg KUg 

Oa ySv* 

where ive have set — — — — = n = constant 

Ga a/2" 

We shall introdtice, In analogy to what was done 
iCarinan for the flat plate, a new varialile z: 

z = ■ §- 



whence 

Y* = 



and 5 = § -5-^ (22) 

Denoting the constant ntunerical values entering into 
the expression for the integrals (21) "by letters: 

5 , '14 _ 23 601 , ^.4. ' 

6 " ^' = ^' ^0 " ^ 1800 = ^' ^^^^'-^^ 
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O 
0 

,o6 



'^^^ = l™(".-^)=^^ (...) 



J 2E^ Ug^ V ^ z / 2K: Tig 

o 

J 2i[^ ug2 ^ S 2 2 > 2^ 



o 



The derlyatives of these integrals wxtii respect to s will 
"be : 

-kf-'^ -I If - A . I) 

G 

08 

_d_ / ,p u du 

a-s 

o 

rlydy=^~^; -i, e^M2--azH^^; 3- e^^ ( 21-.^ 2a- 2)z-^ 

9sJ 2K Ug^ ds 2K ds Ug 

o 

I- /"u^ydy = 1-5 e^^ [z^ - (2a - l) z + ^ ^ a] 

o^,j ii as 

0 

On the right-»hand side of these equations we may write the 
total for the partial derivatives since all the variahles 
depend only on s. 

Suhsti tuting all the values ohtained into expression 
(-18), dividing and collecting similar terms, we find: 
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K ^ 2© dB K ds ^2; ds K ^2 > 



.2 



- cos e -^p- 0^ ■ I (2az-:^21)+a) 7-*^— 7-7" + 



ds u§ ds 



Id cos e . 

4. . ^ ^-.^ . (az-b) 

cos 6 ds 



r2„ 3 



-r. r 



(23) 



The terms containing cos 6 were grouped together so 
as to ena'ble the effect of the term y cos 6 to "be more 
advantageously studied. 

r 22^1 

Let us multiDly the equation by the magnitude | ^ '^rT 

and write the derivatives of Ug , f , and cos $. in the 
form: 

1 dUg _^ d In Ug 1 dr d In r^ 1 d cos ,6 d In cos 6 

Ug ds ds ' r ds ds ' cos 6 ds 

We then obtain: 

ds ds ds 



4. R^^S^BA^^ ^^z 
2K uo r ® 



.dz d In Ug . ^.d In cos6l 

(2az-2'b+a)— +az — 7 ^+ (az-h 

ds ds ds J 



^ Ug (24 ) 



The above expression also serves as the fundamental differ-^ 
ential equation for finding the relation between 25 and 
s, and thus also 6 = 8(s) and Tq = Tq(s). 

It is useful to introduce in the equation in place of 
Ug , its value obtained through the pressure distribution 

on the hull, which value may easily be found experimental- 
ly* This is the simplest method for finding the relation 

Ug = Ug(s). 

Following K. Fediaevsky (reference 14), let us set 

2 

(^^^ = f » where V is the forward velocity of the body 
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or the velocity of the flow' at infinity. 

from the Bernoulli equation p + E^^^ = td + £!„., 

" " 2 *o 2 

where is the pressure at a great ^dl stance from the 

l3od.y, we ohtain after dividing by £™- 

The variation of 7 — # along the hull is easily deter- 

t P V 

mined e xp e r i m e n t al 1 y . 

Making the change indicated above, we find:: 
Uo = Y In ug = In V + In ^/f ; d In Ug = d In 

Suhr/ti tuting in equation (24) and introducing the 
Reynolds Kumber 

YL YL 
Re = ~ = 



^ n C2 ^ 
whence 



YL 

n .= 



Ee C3 ^/2 
we obtain: 

L cos e z^e^^ f /o '^^^ ^ ^ ^ ^-^ ^ . 



Cq v^2 r ^/f 

. , \ d In CO s 6 
(as-b) — 



ds ds 



ds 



Re Cs ^/2 r- r s 

— — Jf (26) 



This differential equation of the first order with re-^ 
spect to 7, han been set up for a figure of revolution* 
For the more simple case, i,G., th:it of the infinite cyl-^ 
indrica,! wing the equation is considerably simplified since 
r then becomes infinite and the last two terms on the 
right-ha,iid side vaiiish. Thus, for a wing, we obtain: 
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as 



e^CAz^-Bz+B) ~ + e^(2A2^>-Bs) — -~ =s ^.-^—-.^ _(27) 

as as L 

where 1 is the wing chord. 

For the even more simple case of the flat plate situ- _ 
ated in a flow at zero angle of attack, the second term o3a" 
the left-hand side of expression (26 )• likewise drops out 
since there is no pressure gradient along the plate 



(-■■-^-'^ = o\. The velocity just outside the boundary 
layer Ug = V constant, that is, f = 1> 

Thus, for a plate, we obtain: 



e^CAz^ - Bz ^ B) ^ = (28) 

dr L 

where x is the variable distance from the leading edge 
of the plate. Introducing in the above ■ expression the 
Reynolds Uumbor R^c computed for any distance x from 
the leading edge, i.e. : 

we obtain: 

dz ^ Re Cs /2 
e (Az- - Bz + B) -zr- ^ — — . 

or 

e2(Az2 - Bz + 3) dz = Oa Tid Rx 

Integrating the- above equation between the limits 
Rq , R^ , and Zq , z, we obtain: 

1 

Rj. ^ Eg =; — - / (Az^ - Bz + B) dz 

Ca -Jl J 

Zo 

wliicli gires, after subst i tuit ing the numerical Tallies of, 
the coefficients A and B and talcing the coefficient A 
outside the integral sign: 



5e-°' 

J^X *^ S-n = ^' 



^ f'(,s^ m , ^ ^) dz 
/ V 15 15 V 



6K^ ^2 Jz^ ^ 

^It should be noted that, in contrast to the equation for 
the figure of revolution, the wing equation may be applied 
for any angles of attack, the latter affecting only the 
character of the function = f(s). 



24 E.A.C.A. Technical Memorandum Ho. 842 



i*e,, we have olDtained the equation for the flat plate in 
the foria given hy larman (reference 4)* 

3. IJITIGEATIOH OF THE EQUATIONS OBTAINED 



It Ib eyldent from even a superficial glance at the 
equations for the figure of revolution and for the wing, 
respectively, that they cannot he integrated hy any known 
method. We have seen that integration is possible only 
for the flat plate. We shall therefore resort to a graph*^ 
ical integration using the method of success^ive approxi*-* 
mat ions* 

W e s e t : 

e^(Az2-Bz+3) ^ cp^; e^(2Az^-Bz) = cp^ ; e^(Az^-Bz) = ^3 

and modify somewhat the expression in hrackets in the last 
terms on the left-^hand side; 

,^ \ ^-z d In y? > , s d In cos 0 

(x:a2-2h4-a) ~ + az -^-^^ — + (az-b) - — ^^^^^^^ 
ds ds ds 

whence^ after grouping the terms, we ohtain; 

, _ . „ - . ^ , . ^ d In cos 6 

a z « ^ . 



\ ds ds ds / ds 



ds 



or, writing the sum of the derivatives a:S the derivative 
of a sum: 

az ™ [2z+ln( //f cos 

e)] ^ A r2(2h-a) + h In cos 8] 
d0 ds 

Suhatituting the a^hove in our expression (25), we olDtaln: 

cp i:^ + cp J:_jLJ2_.t^ + op d I n r ^ L c os 6 

^ ^ '^s ^ ds • 2S& K Ca Vf r 

/az3 X [2z+ln( ^/f cosG)] -z^e^^ ~~ [z(2h-a) + -b In cos 6] 
i. ds ^ ds 

Re Iv^ Oa ^ 



Multiplying the equation "by ds , integrating "between the 
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limits and s, stnd setting 



Re e"" Oa ^ ^ 

k = — — — — ~ ^— =: constant 

L 



we oTDtain: 



z g . Iny f In r 

j d2 =. k ds - 1 d In Vf - p cp^ d In r 

^0 ^o Ixijlo Xti-VQ 

2.z+ln(y" f co s Go ) 
2 P 

^ al„ / Gos_^e. ^2322 d [2 2+ln(/f cos 6)] + 
(22Q-Mn( cos e) 



2(213-30+1:) In cos 8 

5r / ^^-^ z^e^^d rs(2t-a)+l3 In cos 6] 
z Q ( 2 Id - a ) +1) I n cos G ^ 



(29) 



where the letters with suhscript o denote the initial 
values of the yariahles. 

?Phe integration should rigorously start from the point 

Of transition of the laminar into the ttirhulent layer, i.e.»^ 

in accordance with the notation of K. K« IFediaevsky (refer-- 

ence 13) from Sq == t. The position of this point may he 

determined using the so-called critical Heynolds Numher 

Rg ^ which generally lies within the range: 
ic r 



1,600 < Hfi, _ < 10,000 



Using this critical Reynolds lumher Eg 



kr " ^ 



and computing the laminar "boundary layer, starting from the 
nose J hy the laws determined hy the two-term formula for 
the Telocity distrihution (reference 13), it is easy to de- 
termine Sq = t, which is the lower limit* of integration 
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of expression (29). If the hypothesis of conserTatlon of 
momentxim at the transition point is assumed, we should oh-* 
tain: 

In view of the fact, hpweTer, that oscillates within 

a very wide range and depends entirely on the degree of tur- 
hulenc.e of the flow, the ahove operation for finding Sq 
can "be carried out only with great difficulty. 

On the other hand, as tests have shown, to neglect the 
laininar portion - i.e., to assume that the turhulent layer 
"begins at the nose itself, has a negligihle effect on the 
fina,! result and this is particularly true for full-^scale 
hodies, for which the relative length of the laminar por-- 
tion is negligihly small* 

Assuming therefore tha>t Sq = 0 , we shall have at 
the nose 6=0 and Ur. = 7 ^^'Z = 0; whence froin (22): 

0 Q 0 

^o ^""-f - 0 

"that is^ Zq = 0. Wo likewise find: 
Jf^^Q\ lnv^o=^-oo; rQ=:0; In r^^-cx^; fjl:^ cos£o):=-a> 

The initial value cos Gq depends on the shape of the nose 

of th^ hull. fhe tangent to the generating line of the 
figure of revolution is generally perpendicular to its ax- 
is, i»e#, cos 6q = 0* The integral on the left--hand side 
of eouation* (29 ) may he integrated- Integrating hetween 
the limits 2 = 0 to 2: :± z, we obtain ^ dohotihg the in-" 
tegral hy cp4 , 



cp^ == J cp^ dz = j ik%^ ^ Bz + B) e^ dz = 
o ^0 

^ [Az^ - ( 2A+B.) z+2 ( A+B )] -.2 ( A+B ) ( 30 ) 

Substituting the initial values and the integral obtained 
into our equation (29), we find 



H.Ji.G.A, Technical Memoranaum ITc; 842 



27 




In T 




= IC 



cpg d In A/f 



cpjj d In r 



o 




The equation written in the above form is very suitable for 
performing a numerical integration "by the method of suc- 
cessive approximations. The integration procedure \Ye have 
proposed and successfully applied in practice will now he 
described* 

First of all we consider the "zero approxlBatlon*' 
2 ~ f(s). This zero approximation may be obtained by con-* 
structing 6 = f(s) from the equation of K. E. Pediaevsky 
and solving the transcendental equation: 



Having the curve z = f(s), the zero appro3cimation ^ 
we plot the functions under the integral signs on the right- 
hand side against the corresponding expressions under the 
differential sign and perform the graphical integration up 
to the points corresponding to the chosen values of s and 
thus compute the value of cp4 from formula (31), With 
these values of cp^^ we find the next approximation for z 

and ropeat the entire operation the process being contin- 
ued until values of z are obtained that differ very little 
from the preceding values. The latter values of z give 
the integral curve, G-enerally n n agreement to the second 
decimal place (which is entirely sufficient) is found af-* 
tor repeating the process three to four times. 

The integration may also be performed not with respect 
to s but with respect to x, for which purpose there 
should then be substituted in expression (31), ds dx/cosG* 
The above integrating procedure has the advantage in that 
graphical differentiation is entirely dispensed with (ex- 
cept for finding cos 6 f{s) in the case where the shape 



z e 



n L 
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of thp generating line of the hull is not given analytic-* 
ally s,nd where the possihle errors have a very small effect 
on the final result due to the smallness of the last two 
terms in comparison with the preceding terms). This cannot 
"be said of the method proposed hy Moore who, neglecting 
the term y cos G in the expression r* = r + y cos 6 (this 
cprrespoxids to the vanishing of the last two terms of eq.ua** 
tion (31)) succeeds in comhinlng the integrals 




-~ 00 —do 



into one hy introducing a certain differential operator 
a,nd making it necessary to perform gra,phica,l dif f erent ia** 
tion\ In view of the fact that graphical differentiation 
is much more difficult than integration, we "believe our- 
selves justified in sa.ying that our method is superior to 
that of Moore. 

The actual computation was simplified hy: 1) deter- 
mining the functions , CP3 , etc., not from tahles hut 

by making use of functional scales of such length that z 
could "be read, off to an accura^cy of three decimal places 
(.these scales are given in the appendix); 2) performing 
the graphical integration with the aid of the integraph 
manufactured hy the Coradi firm of Zurich, and which at 
once permitted the de termina,tion of the ordinates of the 
integra-1 curve. The actual computatipn was thus reduced 
t o a mininuip. 

The computations carried out using the complete eq^ua-* 
tion (31) showed that the terms ln(«/f cos S )' and h In 
cos 6 under the differential signs in the last two terms 
are. very small "by comparison with the terms 2z and 
(2d - a)z so that they may he freely discarded.^) This is 
particularly true at large Eeynolds Immhers Re, in which 
case z strongly increases. We may then comhine the last 
two integrals into one: 

^ ^„„. ^ _ „„: 

This is true except near the nose and the tail where 
In ^1 cos G approaches - qp which fact, however, only 
slightly affects the results. 
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2z+lii( -/f cos 6) 
a 

/ cosi_. 2 3 e^z d [2z+ln(yf cos 6) + 

_ 00 



z(2"b-a)+T3 In cos 6 

El cosj, zae^^drzCPlo-a)^-]!) In cois G}®- 

^ 03 



aii. 




2z^e^^dz 



2k 



COS _G_ 
r ^ 



z^e^^ 



( 



- 1^ 



l\aK^ cos 6; 

0 



Here, too, as was clone pr e viouslrr , we introduce for convea 
ience in cornpxitat ion , the function 0^ ^ [^^^ '^^X^ " 
under the diiferontial eign. We thus olDtain; 



, / ^ 2kr ^/f ^ 



,2 s 



2 S *^ Z 



(f - ^)j 



dz = 



-e 



2 2 



s.3-(l+jYfz2-z+^) 



a cos £ 
2kr ^/f 
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where xre have set 

5 



and the initial value: 

The final equation for the determination of z = f (s) will 
thus loe: 

s In ^/f In r cp 
cp^rrk / /f ds - / cp d In Vf - / cp d In r - / — T" dcp 



" " ^ , / ' J 2 k r 7f 



5 



(S2) 

fhe above expression is the fundamental working formula in 
the most eonvenient form. 

Having the relation z = i(s) the remaining functions 
of interest to us may readily "be found. Thus, .according to 
(22), the thickness of the houndary layer is expressed "by 

6 = a „ii„^- lif (33) 

K ug S©.. ^2 K 

The shearing stress di s trihution is found with the aid of 
the expression 



K Ug 



whence 



= t: (34) 



The ratio of the shearing stress to the dynamic impact at 
infinity is: 



4 p V 



The coefficient of total frlctional drag for the entire 
hull •DExsed on its voluiae is 
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\?herG X" is the resultant of the fric.tional. forces, ancL 
U, the volume of the hody, is found hy integrating the 
shearing stress distrihution. Wo have: 



St 

X" = / To COS 6 2Tr r ds ^ 2tt p K 7' 



.7 



rf cos 6 



ds (36) 



where s-^ is the length of the generating line of the fig-* 
ure of revolution and thus is given hy 



°u = 




rf cos 6 



--- ds' 



(37) 



We preferred to inte.g;rate with respect to s rather than 
with respect to x as being more convenient since vre have 
the relation z - f (s). If we had performed the integra'^ 
tlon with respect to x from the very "beginning, then it 
would have heen convenient to substitute also in expression 
(37) cos e ds s dx. 

lii concluding this first part, we may ohserve that if 
wo had made use of the MkuradsB-Prandtl formula (6") in 
place of the Zarman formula (12) as the fundamental velocity 
distrihution law, the computations all would have come out 
tho same except that in place of the coefficients: 

.A, 5/6; .B 14/9; a ^' 23/60; h = 6OI/I.8OO 

w G s ho ul d have ohtained: 



A^ r: X; B« 2; a' - 4; h* = | 
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PAET II 

1. COMPtJTATIONS FEOM THii EQUATION DIRITEI) IB PABT I 

Wb decided to apply the computations froto the equation 
derived In Pa,rt I to the hull of the airship "Akron" so as 
to compare the results with those oTDtained from the excel-* 
lent tests for the determination of the boundary-layer 
thickness and shearing stress dlstrihut ion conducted "by 
Freeman (reference 6) on a l/40-scale model at a mean 
Reynolds NumlDer Re = 15.88 X 10^, 

The computation program included the determination of . 
the houndary-layer thickness> the shearing stress or skin 
friction per unit area distrihution^ and the coefficient 
of frictional drag . at four Reynolds Ifumhers: 

15.88 X 10®, 79.4 X 10®, 251.0 X 10®, and 684 x 10® 

The first number corresponds to the tests of Freeman, and 
the others were selected for the purpose of making an easy 
comparison of the results obtained with those given by 

K« Fediaevsky (reference 14) who, using the power law 
for the velocity distribution in the boundary layer, as-* 
sumed that to the above Reynolds Humbers there correspond*^ 
ed the exponents l/?^ l/S, l/9, and l/lO in the for- 
mula for the velocity distribution, and compared his re- 
sults with those given by the tests of Freeman. Thus, 
considering the logarithmic law as more nearly expressing 
the actiial velocity distribution than the power law, it 
was possible indirectly to check the correctness of the ref- 
lation n = f (Re) given by Fediaevsky. 

The values of the universal constant K and the num- 
ber Gg were taken as K = 0»392 = ? . 375 which cor- 
respond to the form for the local friction formula 




3.6 + 4.15 log (Rg Jo'f) 



The constants in the above formula are those of Karman and 
are in very good agreement with experiment. 

The relation z = f(s) was found by solving the funda- 
mental equation in the form (32) by the method of success-* 
ive approximations, the data of K» K. Fediaevsky being as- 
sumed as the first approximation. The operation was con- 
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tinued until the difference in the values of 2 giYen "by 
two* successive approximations differed only in the third 
decimal place, this corresponding to an accuracy in the 
determination of z of 0*08 to 0.13 percent. This degree 
of accuracy, which may appear as too greet, resulted from 
the fact that, in the first place, we decided to use great 
accuracy for the first calculations so as to have confi- 
dence in the results; and secondly, we wished to determine 
the 'boundary-layer thickness 6 with the greatest possihle 
precision-* The values of .6 change very sharply with a 
slight change in z since they depend on the very "sensi- 
tive" function ze^ (33). Thus* the maximum possihle er-^ 
rors in the determination of z give errors in 5 of 1,3 
to l.l percent* {Generally the required accuracy in the 
computation of z was obtained after four or five apprbx-* 
imations. The detailed procedure in carrying out the com-^ 
putations (for practical application) is given in a sup-* 
plement . 

2, DISCLOSURE OF DISOaEPAHOIlS IN THE RESULTS OBTAMED 
FOR THE THICKNESS OF THE BOUNDAET lAYER 



Halving obtained the relation z = f:(s), we construct** 
ed a diagram 6 = f(s), figure 2 (topmost curve). Expect-* 
ing to ohtain with the logarithmic law much oetter agree-* 
ment with experiment than with the power law, we were very 
much surprised to ohtain large discrepancies. For the en^ 
tire length of hull the values of 6 computed hy our pres-^ 
ent theory came out approximately 50 percent greater than 
the eorresponding value computed hy the seventh-power law 
and the test data, the values of 6 from the se venth-^power 
law giving good agreement with the test results as stated 
hy Freeman and as shown in figure 2. 

Our greatest surprise was the very large deviation of 
results from those obtained on the "basis of the seventh- 
power law. It would have "been natural to consider what 
the computation results would "be upon application of the 
velocity-distri'bution law of Prandtl-Nikuradse : 

— In — 
K y 



Since, however ^ using the above formula would have in-* 
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volved a large amount of computation work and taken up too 
mucli time, we decided to make the comparison of the compu- 
tation results for 6 - f(s) from the Earman--*Mikuradse 
formula and the seventh-f.power law, respect iTely ^ for a flat 
plate for i?hich the integrals may "be directly eraluated. 
There \70.s every reaso^i to expect that "feb© relation hetween. 
the positions of the curye^ 8 = f (x) for the flat plate, 
according to the different lairs, would he of the same char*^ 
ac'ter as ohtained for the figure of revolution^ For the 
flat plate the data chosen were: length' L = 5^39 m, ve- 
locity ¥ = 38.45 m/ s , corresponding to a Reynolds Number 
Re = 15 •88 m/s that is, data approaching those of the 
"Akron". The computations were conducted, according to the 
formulas obtained from the fundamental equation (26) after 
i nte.gr at ing, setting r =; co and f = 1. 

■ With the ' veloeity-distrihution latr of Karman , 

cp^dz^cp^^ez [A2^-(2A+B)z+2(A+B)] -2{A+B) = Re x 

0 . (38) 

and i7ith that of Mkuradse-Prandtl : 

r /p 

/ <p\dz=(p%=e2 [z2--4z+6] - 5 = Ee x (39) 

o 

L „ 

whence 6 === — — ze^ 

Ee K V2 

and from the known formula for 8, according to the 'sev-^ 
enth-^power law (reference 13), 

1/5 ^1/54/5 X 

5 - 0.37L Be x (40) 

The values of and K were taken to he the same as 

for the "Akron". 

The results are shown on figure 3* As may he seen, 
"both logarithmic laws give values of 8 much larger than 
are given hy the seventh-power law, the likuradse-Prandt 1 
formula giving results nearer to the seventh-*power law. 
The attempt to vary the value of 0^ (Moore takes Gg = 
10,84) gave no particula,r results as may he seen from the 
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f igure# It was thus possible to conclude from the ahove 
that the reason fpr the large di sagreeinent in the results 
of the present theory with^the se-ven th-power law and hence 
also with the test results^O was to he looked for in the 
theory itself. Very recently we roceived the work of 
Moore (reference 7) where, In his figure 2, which we re^ 
produce in our figure 4, we see a further conf irmatioli of 
our conclusion. Here the results of the (S) mputation of 

6 for the "H.P^L^ long model" at == == 10® (where 

U is the volume of the model) according to the theory in 
which the Kikuradse-Prandtl formula (5") is applied (de- 
noted hy ) give the same order of discrepancies from 

the seyenth-power law (denoted "by 67) as in our case. 

The results of Moore (fig. 5) oix the "Akron" where we see 

a much "better agreement with the tests and the seventh-' 
power law than was ohtained hy us, therefore appear all ■ 
the more strange. Since the identical theory was applied 
to the two models, there should have heen ohtained a dls-» 
crepancy of the same order and this evidently was not so, 
as may he seen on comparing figures 4 and 5. 



3. EXPLANATIOi OF THE DISOREPANGIES OBTAINED 



Bearing in mind all that was said ahove, we decided 
to consider somewhat more carefully the reason for the un-- 
satisfactory results ohtained hy our theory. Up to the 
present we have considered the agreement of the results 
with the experiment applying the theory of Karman in its 
original form. We shall now consider the velocity distri- 
hution itself. We assumed at first that the velocity dis- 
trihution in the boundary layer is computed hy the loga- 
rithmic law: 



and the local friction hy 

= a + h log (Eg/Of ) 

V Of 



Kernpf (reference 10 ) has shown that for a flat plate the 
seventh-power law agrees very well with test results. 
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We shall now attempt, using the excellent tests of 
Freeman, to check the correctness of these laws- We have 
at our d-isposal the tahles of u/u§ = f(y) for the Tari- 

ous sections of the hull (placed at the end of the report 
of Jreeman) and also the experimental curve of shearing 
stress distrihutipn (the curve 0 = '^o/^o fig« 8 of 

Freeman ^s report), The change-over in the computation 

u - Ug 

from the values u/uii to the values -^—.^-^ easily 

• Ho 

effected thus: v p 



u u 



"^o ^ ni 




where ^ ^ p Y 



^o 



If the formula for the velocity distrihution Is correct, 
then "by drawing the test curves, using as coordinates 

^ and In (l - /l - Z W / 1 Z ^ we should ohtain 




p 

a grouping of the points ahout a straight line passing, 
through the origin of coordinates and having a slope de-^ 
termined "by the maghitude of the universal constant K* 
We pointed out at the heginning that the current method of 
determining K was hy tests on the drag. For this reason 
the method proposed hy us for the determination of K from 
the velocity distrihution appears to us to he original and, 
as we shall, see helow, leads to results of extreme impor- 
tance. 

If we consider further the formula for the local frlc^ 
tional Intensity distrihution as true, then wa should oh-^ 
tain a straight line on constructing the friction curve 
using as coordinates: 



„ Jf 

^0 



7^ 
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an d 



log (Eg /Gf) = log 



and haying a slope that determines the coefficient h and 
heiice also the same universal constant Km 

Figures 6 and 7 ' show the above curves thus plotted. 
The points — — - do not lie on one straight line "but 



arrange themselves along some curve showing that law (6*) 
is in general not strictly true as, of course, is to he 
e3Cpocted for the regions near the wall (at the right of 
the dia.grains ) where the viscosity has not heen taken into 
account and. where the curve very appreciably slopes down— 
ward* In the region, however, sufficiently far removed 
from the wall, it appears possible to draw a straight line 
through the points, which fact shows that in this region 
tho universal law is valid- 

The result of the computation of the universal constant 
from the slope of this straight line led to the value*^) 
K - 0,214 in place o.f the value .0,4 generally assumed* 
Gonsidering now, figure 7, where out of a very small number 
of tost points many dropped out (in the region from 9 to 12 
foot from, tho nose, according to Freeman (x/L from 0.4 to 
0,5), where the curve t'q/^o - ^(^) takes an improbable 
dcwnwvard slope). Five of the points, however, aline them-* 
solves sufficiently well on a single straight line, from 
whose slope and intercept the values a = 3.4,- b = 3»77 
are obtained., in very good agreement \yith the values as-^ 
sumcd by us In the computations, namely; a = 3.6, b 
4.15 and from which tho value K i= 0.389, which is gen-^ 
orally accepted, is obtained* 

Having obtained these results from the tests on the 
"Akron", it was natural. to reconsider the classical tests 
of Hikuro.dse on pipes, in the light of these results those 
tests being the only generally accepted source for furn-* 
ishlng the universal constants. We made use of the data 
of tabic Yll of ETikuradse on. smooth pipes (reference 3). 
The plot is given on figure 8. Contrary to what was ex-^ 
pectcd, In the r egion sufficiently far removed from' the 




The curve corresponding to K - 0.392 is given on the 
diagram. 
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wall and for sufficiently large Reynolds Ivfurn'oers, the same 
smaller value of K, namely, K = 0»3G6, that is, mucii 
less than ivhat is currently assumed, was again detained 
also in thi^ case, Prom the same test data, "but applied 
to finding the cvrag ooefficient, values of Z ivero oi-^ 
tained near 0#4 (the corresponding straight veloci ty-dis-^ 
trihution lines are shown on figp 8). We are thus led to 
the conclusion that the "universal constant" hj no means 
^PP^^^^ to he universal. In all cases of tests on the ve*^ 
locity distrihutiori the constant comes out much less than 
in the tests on the drag, 

Besidfs the tests of Nikuradse, we have also employed 
the data of Hansen (reference 11) on a flat plate. In 
this case an unusually low value of E == 0.175 was oh- 
tained. 

It was also of interest , naturally, to work over the 
tests of Freeman, Hikuradse , and Hansen in terms of In s/y, 
that is, to check the ITikuradse-Prandtl formula. fhe 
curves have a form analogous to those just analyzed, and 
the values of the constant In the region far removed from 
the wall, a.GCording to the tests of Freeman, Hikuradse, 
and Hansen are, respectively, 0.267, 0,350, and. 0»22, 
An excellent explanation of this phenomenon after it was 
disclosed, has heen furnished hy Professor L. S* Loitsansfcy, 
who states that the reason for the considerahle variation 
in K is the result of neglecting the effect of the vis- 
cosity in the theory. Actually, the application of the 
theory of similitude gives a logarithmic law only under 
the conditions where the effect of viscosity is negligihle* 
In those, regions far removed from the wall, where the effect 
of viscosity is negligihle, the points follow the logarith- 
mic law with sufficient accuracy with corresponding actual 
value of K (0,2 0»3). As the wall is approached, the 
effect of the viscosity increases and the points hegin to 
deviate from the straight line in such a manner that the 
derivative 



v* 



decreases continuously (corresponding to a continuous in^ 
crease in Z) and in the regions near the wall the in--- 
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clination.of the curve (as may he soeii on the diagrams) 
corresponds to a yaluo of K z 0.4, 

The latter circumstance is particularly noteworthy* 
Aeeording to the classical theory of Karman, the friction** 
al drag la^ is obtained froia that of the yelocity diBtri-^ 
hut ion by applying the latter to the conditions on the 
■boundary of the laminar sublayer (which process , of cour se ^ 
is formally unjus'tified since the effect of the Tiscosity 
is still very appreciable) - i.e., to the region where the 
value K ^0.4 should be taken. Since the value of K 
has, rp to the present, been determined exclusively from 
tests on drag, the value always obtained was near 0.4, 
i.e., this value' represented conditions at the boundary 
of the laminar sublayer .in the region very near the wall* 
This result was likewise obtained by us in figure 7, 

■ There was thus obtained a pa,radox: The usual method 
of the det ermina.tion of K by studying the variation of 
the resistance -of pipes with Reynolds Number gave values 
of K for a region for which the logarithmic law, strict-^ 
ly speaking, was not applicable and the region where " the 
law was applied - that is, at some distance from the vrall 
was not assigned the corresponding va.luo of E which, as 
we now SCO, is much less than the generally accepted one* 



4* IMPaOYEMElTT OF THE LOGAaiTHMIG TEEOHY 



!From all that ha.s been said above, the following prac-* 
tiCc?vl conclusion may be derived. Since v/e do not at the 
present time possess any velocity~-distrib^^tion law that - 
t-a,kes into account the effect of viscosity, we are con*^ 
strained to use the logarithmic law of Karnan, rendering, it 
more accurate by assuming for the region fa-r removed from 
the wall a value of K - 0.2, and near the wall K 0*4, 
We cannot, however, determine definitely -just whore the 
value of K under?^oes a discontinuity (corresponding to a 
hypo thetica,! limit where the visco'sity ceases to have any 
effect). Q) Moreover, the breaking up of the interval y/s 
from 0 to 1 into two parts \vould introduce .codlications in 
the computation. We therefore decided that, with very 



■It is interesting that, according to the form of the dis-* 
tribution of the points u-ug/v* for the "Akron", a sharp 
discontinuity is observed at the value y/6 == 0.19 (fig. 
6). 
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little error, we maj^-^arrive at the following compromise: 
for the veloclty-distrHsutipn formulas' to eonsider for the 
whole region Kx - 0.214, anii for the resistance formula 
K 1^ 0.392, There ±k no douht that the Yalue Kx is af-* 
fected. also hy the not entirely correct procedure in pass-^ 
ing from the internal to the external prohlem. It is q[uite 
possihle that for other cases of the external prohlein dif-* 
ferent values of Ki would he ohtalned. The only thing 

that may he affirmed with any certainty is, that, in every 
case Kx trill he less than K* This is confirmed in all 
three cases of flow considered here ^ 

. Having assumed that the reason for the very large de-^ 
viatlans in the computed values of 6 from the test val--^ 
ues is the incorrect choice of K in the formula for the 
velocity di strihut ion , we decided to reconsider the theory 
on the basis of ^ the ahove compromise assumption that leads 
to a choice of two values for K» 

Assuming the value Xi = 0.214, we obtain the follow^ 
ing expressions for the integrals 
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i.»e,, vre finally, after stilsstitutl ng 8 from expression 
(22), obtained different cpef f icieijts in pla&e of A, 3, 
a, and "b, as follows: 

IK ^ BE^ aK ^ "bK^ 



The functions cp^^, cpa* CP3, CP4, and cpg in the funda-^ 

mental equation will then haire different coefficients: 
Denoting the new functions hy cp^^^, Q^Si' ^3i' ^4x- 
cpg we obtain : 

cp, =eMAi22-Bijz+Bi)=02 J_ gS^B _K1 ^,^3 _K \ 



(^3 =e2(AizS^B^z)=eZ (2A — z^-B z) 



%^=Q^ Ai z2-(2Ai+Bi )z+2(Ai+Bi ) -2 (Ai+Bi )=e^ 



AK o 



. (2A+B X) z+2 (a+B JLI .2 (a+B ^ 

\ Ki/ El \ El/ El J \ El/ El 



= e2 2! 



(42) 



The final equation, analogous to (32), for the deterinina- 
tion of z = f(s), applying the law with two values of E, 
will he: 



/ ^ii cL2=K^0g -/2 / Jf ds- 



s In 'Jl 

^ ^ - ' q^a d In -y/f - 



00 



In r 



cp 



(43) 



cpj^^ d In 



La 



cos . 9 J 

2Eie2 -/2 Re J r Vf 
%i(z^o) 
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e qual to 



CP . ^ ^ I ^ l) 



Since in the above equation the functions cp. enter 
with different coefficients (in accordance with (42)), It 
would appear necessary in solving it, for us to recompute 
all our functional scales of which we- made use for rapid- 
ity of calculation. This, in the first place, would have 
very much delayed the completion of the work, hut what is 
more important, would not have "been justified since the 
value ICi = 0.214 has not yet "been estahlished as a unl-^ 
versal const8.nt, and on changing it the recomputed scales 
cpg^, ^3^> ^4 3_' ^5 3^ would have hecome useless. 

¥e therefore decided to reduce equation (43) to a 
form in which it may he used with the original functions cpg , 
cp3 cp^, and cp^ as may he done, as we shall see, hy in** 
troducing in it one extra term. 

"We find the expressions cp-^^, ^2 ^ ^ ^Zx' ^4i' ^^^^^ 
cp^. , cpg , CP3 , cp^. From (42) and (28) we ohtain: 




(A:2^-Bz4.B)+Az^0^(^|- -1) 




and similarly 



for the remaining functions: 



^ J:^ L ^ 2A 2^ e^ - l)] 




The function 



cp4^ we define as 




that is : 



o 
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d2 





e 



2 z2 ^2 



To olitain a suita"ble expression for the grapiiical Integra-- 
tion, it is "better not to compute the integral on the right-* 
haiid side of the ahove expression. With regard to the last 
memher of the equation which ta,kes account of the effect of 
the term y cos 6 we decided, instead of replacing cps^ 

^Psi to. compute the scale for cp^ directly (setting 

Ki 0.214)* At some other value of K (for a different 
case) this would have a negligihle change in the final re--* 
suits since, as we shall see "below, the effect of the last 
term of equation (43) is in general very small. 

Let us now suhstitute all the expressions obtained in 
(43) and split the integrals with respect to In ^/f and 
In r into two* We then ohtain: 



z 




e ^ d z 




1 = K ^Ca v/2 7 Vf ds - 



o 



In Vf 



I n 




In r 



In r 




cpgd In r 




e^ z^ d In r 




44 



II.A.C.A* Technical Memorandum lo, 842 



Multiplying the entire expression "by Zi^/Z^, rearranging 
the terms .and com^binlng the integrals of the functions 
e^z^ into a single integral with a complex differential, 
we find: 

in Jf In r 




I B ^^^^ 



ZK^ Cs 72 Re J T ^'^ 

tp5i(z = o) 

We have thus ohtuined an expression that differs from 
(32) hy one additional integral and has different constants 
hefore the first and la.st terms on the right-hand side. 
The presence of still another integral prolongs the compu-^ 
tation somewhat ("by a^hout l/5) hut on the other hand the 
new expression, using the same scales of the functions qp 
permits computation with ©.ny value of ICi* 

The a-hove method of introducing into the equation an 
extra term that takes into account the variation in the 
coefficients cp is of advantage not only when the change 
is "brought about hy the equality of K and IL^, hut in 

other cases as well. In particular, as has heen pointed 
out in Part I, the coefficients may vary as a result of 
the application of a different logarithmic formula for the 
velocity dis trihut ion* Let. us put 

ill 



Prooeeding with the functions cp as ahove, we ohtain fi- 
nally? 
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2+ In rf 

+ A (l - ~|-) ©5!^ 4(2+ la yf ) « 



5x 



2K^ P Cs /2Ee r 7 f 

T/here the change of coefficients in cp'^ may likewise be 
neglected. ^ 

The aoove expression is suitable for computation ^.yith 
our scales using any coefficients A and B. In particu- 
lar, i'f it is desired to use the veloci ty-di stritjut ion 
formula: 

u - U ii ■ 1 ^ S 

^* Ki y 

we mist take: 

a =: 1 : I = 1.2 p = 2 : 1.286 



:-5. REPEATED COMPU TATlOI^fS IT.SIN& BQUATIOH (44) 



'The computations from equation (44) were carried out 
for the same Reynolds Numhers: 

15.88 X 10^ 79.4 x 10^ 251 X 10^ 684 x 10® 

with ar).d without the last term tha.t takes into account 
y cos 6, The method of computation was the same comhlned 
graphic and analytical one, successive approximations in 
the Talues of' z helng taken until a difference in the 
third decimal place only was obtained. 

On figure 9 are shown plotted the changed values of 
z for different Reynolds llumhers against the ratio x/L 
along the "Akton" hull. As may he seen, the effect of 
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taking account of the last term with cos 6 shows up, as 
also may l)e expeGted, only at the tail itself where the er- 
ror in rejecting the term y cos 6 in the expression 

r f - r + y CO s 6 
is a maximum as iaas also heen pointed out hy K« Jediaevsky. 

The strange "bends in the curTes z =: f (x/L) on figure 
9 in the region of the stern result from the fact that the 
function 2 In rf , wi th respect to which the integra/t ion 
in the next to the last term was performed has a wavelike 
form in the region of the stern. 

Having ohtained z == f(s), we plotted on figure 2 the 
changed- thicknesses of the "boundary layer 6 == f(s) for 
Re - 15.88 X 10^. As may he seen, excellent agreement was 
ohtainod with tlie test results much "better than was given 
hy the soventh-power law. Particularly noteworthy is the 
good agreement in the region of the stern (s = 4 5 m) 
due in large part to the fact that y cos 6 was taken into 
account as may he seen from a comparison of the curves 
8 - f (s) with and without the corresponding term account-^ 
ed for. 

A comparison of the shearing stress computed from the 
formula : 

'^o 2K^ f 




according to the original solution Z = 0»392, the "com- 
promise" solution with K = 0.392 and Kx = 0.205, and 

the power law, with the results of the Freeman tests, is 
given on figure 10. 

We see that our compromise solution gives the "best 
agreement with experiment. In -the center portion^ without 
considering the douhtful hend in Freeman ^s curve, our curve 
almost coincides with the test curve. In the tail region it 
may again he ohserved that taking into account y cos 6 
gives a hetter agreement with experiment. As regards the 
very great discrepancy in the nose region, this phenomenon 
is unavoidahle. In the first place, due to the fact that 
the hull contour curvature, \7hich is considerahle near 
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the hose, 9^ was neglected, the theory is not entirely true, 
and in the second place (and this is the more important 
reason) the .accuracy of the experimental det erminat ipn of 
the skin friction, due to the smallness of the houndary 
layer, is not very large. 

From the results ohtained ahove, we couid plot the 

drag coefficient against the Reynolds Uumher. The value 
of the coefficient for each Reynolds llumher is ohtained hy 
coinputing the expression: 

1 

^ / £| a ^ 

o" 

The results are plotted on figure 11, which shows the test 
curves of different models of the "Akron" and also the com-^ 
puted ciirves of Moore and K. K. Fedlaevsky • If^) 

We are justified in saying that oue curve more closely 
represents actual conditions than the other theoretica-1 
curves since it agrees better with the results of the most 
reliahle tests conducted on a metal mod,el in the variahle 
density tu.nnel. The forn drag, i.e., the integral of the 
excess aerodynamic pressure along the hull was not taken 
into account since, according to repeated calculations 
(reference 8), this drag component appears to he very near 
zero* We may thus consider the coefficient 0-^ as heing" 
that. of the total drag. 

Figure 11 also shows two points computed hy our theory 
with K - El = 0.392. They lie considerahly "below the , 

curve with K = 0*392 and Ki == 0.214 as should he ex- 
pected. The curve of Moore lies hetween the latter curve 
and the points. This apparently shows that the velocity 
di.strihution of Nikuradse-Prandtl : 

'^0 K y 



^In introducing the function r = f(s), we introduce not 
the curvature of the contour "but the law of variation of ' 
the surface at which the friction is developed. 

■'■^Wc curves aro given as functions of the Reynolds iJumher 

hased on the volume to the 1/3 power : R^ = -^^-^ • 
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foi-* the condition K = ICi - 0.4 corresponds more nearly 

to the truth than the Ear man formula under the same condi*-^ 
tions. This is likewise confirmed hy the fact that the 



value for Ki than Ei = Qe214. Thus Moore, choosing the 

H ikuradse-Prandtl formula and not taking into account any 
difference l^etween K and Ki, committed a smaller error 
than we did in our first computations. The curve Ij^ing 
lowest is that computed "by the method of K. K. Fediaevsky 
with n = var. 

Figure 12 sho^7S a comparison of the variation of the 
relative thickness of the houndary .layor with distance 
along the hull according to the present theory using K ~ 
0.392 and ICi >= 0,214, with the computations of 

K, !C» Fediaevsky. Although . the curves agree sufficiently- 
well for Re = 15 .SB X 10^, the agreement is far less at 
the higher Reynolds HumDcrs, Prom this it nay he conclud- 
ed that the relation n = f(Re) proposed 'by K • K, 
fediaevsky, does not entirely correspond to an approxlma-^ 
tion to the logarithmic law. We found 'by integrating with 
respect to the curves 8 - f(n) and averaging the results 
for different sections, a new relation n = f (Re) (fig. 13), 
which should more nearly "bring the results of the computa- 
tions of 6 hy the power law closer to the analogous com-^ 
putations hy the logarithmic law. 

in conclusion , we may note that with very considera-» 
hie accuracy, for more rapid computation of the relation 



plied, computing the houndary-layer thickness hy the power 
law- (using the corrected curve n .= f(Re)) and the shea,r- 
ing- stress distribution and drag coefficient hy the loga*- 
rithmic law. !or this purpose, having ohtained the values 
of 6, the values of z are next computed hy solving 

e qua t i on (22) for 2 e ^ : 



u 



plot of the points 




against In — gives a larger 




a "mixed" method may he ap-^ 



ze^ = 



Re E Gs 



8 Jf 



(45) 



after which Tq/^o ^"-^^ ^u computed hy formulas (35) 

and (37). This method of computation gave curve = f (Re) 
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lying only slightly a*bove the more accurate curve. For 
practical purposes this method appears entirely accept-* 
ahle # 

6. GOHGLGSIONS 



The results of our iaayestigat iou lead to the follo^T"- 
ing conclusions: 

a)' The logarlthinic law of velocity di strihutlon of 
Karman, applied to a figure of revolution and to a flat 
plate with the generally accepted value of K 0.4, gives 
too sharp an increase in the boundary-layer thickness arid 
too small values for the drag. 

h ) This Is due to the fact that a ouch smaller value 
of the constant K (Zi ^ 0.2 ^ 0.3) corresponds to that 

region where the logarithmic velocity-- di strihut ion law is 
actually applicahle (i.e.>, where the effect of vdscosity 
i s nogligihle ) • 

c) If the logarithmie law Is applied to the region 
near the wall where the effect of the viscosity shows up 
strongly and where this law, strictly speaking, is not ap-^ 
plicahle, an error in principle is made which it is possi-^ 
hie to compensate hy taking K iz 0.4. 

d) Since a variable K cannot he introduced, it is 
proposed that in the formula for the velocity distrihu*^ 
tion a value ICx « 0.2 - 0.3 he taken, and in the formula 
for the drag, Z z 0.4. 

e) With the ahove values assumed for excellent 

agreement is ohtained with experiment in all respects 
(thickness of the layer, shear stress distrihution , total 

f) Further light on the significance of K for the 
velocity- distrihution law can he ohtained only hy more ac-* 
curate tests on plates, wings, and figures of revolution, 

g) The method proposed hy us of employing two values 
of K for the computation hy the logarithmic law Is of 
course a temporary e:v:pediont to he applied until such time 
as a more accurate law that takes into account the 'effect 
of viscosity - i.e., one correct for all regions of the 
flov/, is substituted for the logarithmic law. Obtaining 
such a law is one of the most urgent of present-day prob- 
lems • 
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APPENDIX 



Some Practlc<al Suggestions viith Regard to 
the Computatipn Procedure 



Por carrying out the computatipn using the complete 
equation (44), we may here present some practieal sugges- 
tions j gained from our experience, to those trho wish to 
make ivse of our method. 

( ) Shoioe 0 f _in ijfci al_e 0 n^ •-•It is "best to he- 

gin the integration from the end of the laminary portion 
of the "boundary layer. Por determining the transition 
point from the laminary to the turhulent layer, it is nec-* 
essary to make use of the critical Reynolds Humher ("based 

on tne houndary-layer thickness) Rg = — ^ and which lies 
within the range 1,600 < Eg < 10,000« Prom the corre- 

spending value of Rg , ^lam ^ ""ug determined; whence 

a.ceording to the transition hyoothesis of Karman (reference 
1=5.)' huxl^ = Ham' 

After this, prpceeding to z according to formula 

(45), the integration Is "begun, the Integrals "being taken 

from the section chosen and the corresponding initial value 
of the function cp^ found. 

As may he seen, however, the choice of the critical 
Reynolds lumher is suhject to a certain degree of indefi*^ 
nitonoss# Per this reason it is iDosslhle, with very little 
error, to take the lower limit of integration at the nose 
itself, sotting z = 0 at s - 0; that is, consider the 
turhulont layer to start at the nose itself. This procod-- 
ure will he m.ore nearly justified . the larger the Reynolds 
Humher, Re - YL/v , taken. Here a new dif f icul ty ari se s , 
however. Tilth the initial conditions z = 0 when s = 0, 
the integrals with respect to In Vf e.nd In r must he 
eva-luated from od (since when s = 0, f = 0 and r == 0 ) 
which, of course, cannot he done graphically. However, as 
computations have shown, the value's of these integrals con- 
trihute very little near the nose on account of the '^weak- 
ness" of the logarithmic infinity* Por this reason the 
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lower limits of these integrals may, \-r±th Tery little error^ 
"be taken at any point sufficiently near the nose (for ex*^ 
ample, the point corresponding to In ^/f S - 3.0) neglect- 
ing the integration up to this- point. The integration may 
also Ido hogun from any point sufficiently near the nose, 
having determined the variation of 6 from the nose up to 
the point hy the method of K. Z. Podiaevsky, iTith n = var» 
This is, in fact, the most accurate method of integration 
assuming 2 = 0 xrhen s = 0. 

(h) Ghoice of the first approximation As a first 
approximation, the relation hetween z and s" ohtained 
from the hpu.ndary distribution hy formula (45) according 
to the method of Km K. Fediaevsky, may he taken, as was 
done hy us* It is also possihle to take as the first ap*^ 
proximation the solution of equation (44), neglecting all 
the integrals on the right-hand side except the first, i*e«, 

cp^* = K C2 ^2 J ds 

o 

as T:as done hy Moore, Having obtained the value of cp^ ■ , 
z ' Is found "by using the corresponding scale, 

( 0 ) Ohtaining the suc ceed ing; a-Dproximati on s • Af t e r 
obtaining the relation z == f(s) from some one or other 
first approximcation , vre find the values of the functions 
cp^, cp^ 5 , lolot the functions to he integrated against 

the corresponding expressions under the differential sign, 
X • e • • 



cpa = f (In /f ) (P3 = (In r) 

e2z2 ^ fg(z+in rf) and == % (qpg ) 



(46) 



and integrate the curves graphically between the limits 
chosen on the meridian line s# Having obtained the val^ 
uos of the integrals, xre substitute them in expression 
(44), determine a netr value cp^" and a new value 2" and 

repeat the operation until two successive values of z 
differ by as little as required for the accuracy desired. 

In performing the integrat ion outright along the en-* 
tiro length of the meridian, however, the following dif-^ 
ficultios are met with. The functiops In and In r 
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generally vary sharply near the no se and tail { » 

r— >0), T.7hereas in the center portion the change is very 
snail, The cp functions, on the contrary, vary much more, 
sharply at the center portion than near the nose. lor 
this reason it is quite impossihle to choo se for the on^ 
tire length of hull the same scales for plotting (46) that 
give inore or less the same acciiracy in the determination 
of • z in the difforent regions. It is therefore suggest-^ 
ed that the distance s l^e "broken up into throe intor-^ 
vals-^^^ and. the integration performed for each of those 
separa.tely» It may appear that this device wotild prolong 
the computation. The computation, hoiTever, will, on the 
contrary, "be shortened ^Thile the accuracy will at the same 
time he increased. The shortening of the s.mount of labor 
will be accomplished in the following manner: 

Having integrated over the first interval and finally 
obtained the real value .at the end of the interval, 

the values z^'^ that were obtained from the first approx*^ 
imation are corrected in the second Interval by multiplying 
them by the ratio z^/z^*)^ that is, by the ratio of the 

true value at the end of the first interval to the value 
of the first approximation at the same point. We thus ob-^ 
tain a.t once for the second interval a better approxima^ 
tion to the actual curve z = f(s). Gomputations ha,ve 
shown that by such a, procedure good values of z were ob** 
tained for the second interval after two or three approx-* 
imations, while for the first interval it was necessary 
to take five or six approximations. The same method is-ap-* 
plied in pa^ssiug from the second to the third interval* 
A better result may also be obtained in passing from one 
interval to another by multiplying the function e^ or 
cp^ corresponding to z of the first approximo.t ion by the 

ratios of the values of these functions at the end of the 
first interval. 

To give an idea of how the plot of (46) looks for the 
complete graphical integrationj we present figures 14, 15, 
and 16^ showing these plots for the first, second, and 
third intervals of the model of the ^*Akron" at a Reynolds 
Number He = 15.88 x 10^ and corresponding to the data of 
the last approximation z = f(s). The points are numbered 



divided the length of the "Akron" into three portions, 
thus: the first from x/L = 0 to xy^L == 0^1 ; the second 
from x/l = 0.1 to x/L == 0.7; and" the third from x/L = 
6.7 to x/L = 1,0, 
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according to the niimlDers of the sections chosen. The 
ta*ble (p«.56) giVes for each of the numhered points the 
distance frorn the nose along the hull and the values of 
f , r , and cos 6. 

The graphical integration of the curTes of a given 
type ^vas performed using the integraph of the Ooradi firm 
of Zurich an apparatus constructed on the principle of 
ATDdank-^Ahakanoyi tch and which hy following along the 
curve, gave the area under the curve "between the initial 
point and any point on the curve. If the ahove apparatus 
is not availahle the usual planimeter may "be used or the 
area found hy summing up the squares under the curve:. 
For ccnvenienca In computing hy our method, we give in 
the appendix the functional scales of cp as functions of 
z . 

Our scales j however, do give not the functions cpg , 
cp3» Cp^ , ' ^"^^ these functions divided T)y K == 0»392; 

that is, 

and 

0,392 0.392 0.392 0.392 

Thus, in using these scales, it is necessary to divide 
the terms of the original equation "by 0.392, 

In conclusion we may observe that, in accordance with 
our computations, the term containing cos 6 may "be. en- 
tirely neglected in the first interval and at the "begin-* 
ning of the second. Only farther along the hull does this 
memlDer show any effect, the latter showing up only in the 
thickness of the houndary layer and in the distribution of 
the shearing stress hut having practically no effect on the 
integrated drag. If the computations therefore a.re made 
merely for the purpose of determining the curve 

f (Re ) , the last term may he entirely neglected. We may 
8.1 so point out that for rapidity of computaticn we may, in 
the fourth term on the right-hand side of the last equa^ 
tion, write.: 

Az® e^ " *Pa ^3 

i.e#, In the form of a difference in the functions used 
in the two preceding terms. 



Translation hy S. Reiss, 
Hational Advisory Committee 
for Aeronautics. 
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Figs. 1,7, 9 
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3.4 3.6 3.8 4.0 4.21^(55/0^) 

Figure 7.- Check of the Zarman 
relation iJL=£ =a-*-'b Ig 

(RgV~Cj) from Freeman's data on the 
airship model "Akron. 




FigTare 9,- Variation of z along the hull of airship 

model "A^cron" for different Reynolds Ifumhers 
(our theory with K = 0.392; == 0.214) , 
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i?i£M.re 2,- Bo-oniry- layer thickness 5 in meters along Imll for the 

l/4:0-scale model of the "Aicron" at He = vL/u = 15.88X10*^. 
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Lofisaritlunie law of 
Karman Ci=7 . 375 ,* K= . 392 
--Ci = 10.84 
jLogeritlimic law of 
' Nikuraise C-^^. i.37.5 
= 10.84 

J Logarithmic law of 
I Karman K= . 3 92 ; K| . 24 
"-Power law n=l/7 



Figure 3,- Boundary- layer thicicness ^lon.r; a flat plate 

according to the various theories, velocity 
of flow V = 38.45 m/sec. 
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WLgors 10.-. Sba&rlng 9tr#s8 dtttrl- rigor© 12»-. Variatioa of relative 

Imtioa along taill of air-- l^otsndaxsr^-layar thickaosa along Imll 
ship model ^ikron^ at Re P of aireMp model '^Akroa** at vario-oe 

1S.88X10^« Beyaolde Htmbere aecordiag to the 

logarithmic aaid erpoaeatlal laws. 
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|Our theory K ^ •392, Ki .214 
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figure II*.. fheoretical aad experimeatal Yariatioa of frictioi^l 
drag of hall of airship model '^Alcroa*^ with B^olds 

Number* 



i.A,C.A. Technical Memorandum No. 842 



I'les.13,14 



1/n 

10 
y 
8 

































































































/ 



























































































































































































































































































































-iTroni data on ^' Akron" 



:: - - K. K . Fedyae vsky 



7 8 9 log Re / 

Figure 13.- Variation of the exponent in the pov/er 

law of velocity distribution with 
leynclds N-um'ber, 




Figure 1*±.-^ Integral curves for interval I. 

Model of airship "Akron" ,Be=l5.88xl06. 



IT,A,C.A, Teclinical Memorandum No, 842 



figs. 15, 16 



o 

X 

M 



5 qpS^/KxlO"^ 




0,6 Inr 



Fi^rare lo..- Curves for interval II, Model of airship 
"A^ron", Be 15.88x10^ 



at) 
\ 

o 

X 



.-12 



t0 
f 



X 



100 



50 





1 


0.2 




10 


'3 


io.4 




L0,5 z+lnzf 
























































30 














































































































-i: 


















































































































p 


















































IS 








t 




JL ugs o 












— 


.2 




) 












2K^C i\/2 Ee 2 Vf 




























17 








































^^ 
































1 








18 






























fa 










ij 


j 




^-^ 




















If 












hi 












n 


-6 


















i /K 






i 










7 










































o- 






1 




































































X 


-,05 0 


.05 


.1 


In V?" 



-i.y -1,7 -1.5 -i';3 -1,1 -/9 
Figure 16,- Integral curves for interval III, Model of 
airship "Akron, 
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